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Abstract 

We present several new asymptotic trace formulas for Jacobi matrices 
whose coefficients satisfy a small deviation condition. Our results extend 
most of the existing trace formulas for Jacobi matrices. 
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1.1 



Jacobi matrices have a wide range of applications in mathematical sciences. 
In physics for instance, they naturally appear in the study of random matrices 
and discrete Schrodinger operators 5 . In statistics, they provide a useful tool to 
study stochastic processes such as the birth-death process and random walks. 
In classical analysis, they play an important role in the study of orthogonal 
polynomials [51 H3]. 

The spectral properties of J(a fc ,b fe ) are well-known and can be found in 
many texts; a good reference is e.g. [6]. If M denotes the quantity 



M = V3 



max \af\+ max 16? 

Ki<k Ki<k 
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then one can easily verify that a(J(a fc , b fc )), the spectrum of J(a fe , h k ), consists 
of k real simple eigenvalues lying inside the interval \—M,M\. 

Except for few special cases, it is in general impossible to obtain explicit 
expressions for the eigenvalues of an arbitrary Jacobi matrix. However, an 
important special case for which the spectrum is explicitly known is the case of 
tridiagonal Toeplitz matrices. These matrices have the form 



T k (a,b) 



(1.2) 



fa b ••• \ 

b a b : 

Ob a '■■ 

: '■■ '•• '•• b 
V ■■• b a J 

for some net and b > 0. The eigenvalues of Tfc(a, b) are well-known and are 
given by 

= a + 26 cos f^ 71 ^ for j = 1, fc. 

Using Riemann sums, it is then straightforward to derive the following asymp- 
totic trace formula 

k 



lim ^-Tt&ceWTda, 6))] = lim i V cj> (a + 26 cos ( J -tx\ \ 

fe->oo k fc->oo k ^ \ \k J J 

i r 

= — I 4>(a + 26 cos x) dx 
it Jo 

for any continuous function 4> on the interval [a — 26, a + 26] . This trace formula 
was successfully used to obtain new ones for Jacobi matrices that are small 
perturbation of a Toeplitz matrix. For instance, this approach can be used to 
derive Nevai's result [11] on the density of the zeros of orthogonal polynomials 
that belongs to the M(a, 6)-class. 

In this paper, we continue our investigation initiated in our recent work [2]. 
We start by deriving two asymptotic trace formulas for the moments that extend 
the ones given in Lemma 2.3 of [3J. In the third section, we derive our main 
results. Namely, if we assume the bounded sequences (a fc )fc, (b fc )fc to satisfy the 
small-deviation conditions of Definition 12.11 and if (a , b fe )fc is //-distributed as 
in the sense of Definition 13.11 then we have 

7T 

4>(x + 2y cos t) dt d(i(x, y) (1-3) 



lim -|-Trace[0(J(a fe ,b fe ))] = ~ I [ 

fe^oo k 7T J[0,1] 2 Jo 



for any <fi E C[— 3, 3]. We also present similar results for unbounded sequences. 
In addition, we show in the fourth section of the paper that for any given 
sequence (a fc )fe, (h k )k satisfying the small-deviation conditions of Definitions 
2.11 and 12.31 one can always find a subsequence and probability measure fi for 
which a trace formula similar to (jl.3p holds. 
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2 Moments 

We start by proving a trace formula for the moments of J(a fc ,b fe ) when the 
sequences (a. k )k and (b fe )fc are bounded and satisfy the small deviation condition 
below. We also present results when the sequences are unbounded. 

The first thing consists of defining in precise terms what we mean for a 
sequence (a fe )fe with a fc = (af , a^) £ K fe to satisfy the required small deviation 
condition. 

Definition 2.1. We say that the sequence (sL k )k belongs to S if it satisfies the 
following two conditions: 

ft) EiJi 1 ki+i - al t\ = °( fc ) as k -> °°> 

(ii) < a k < 1 for all i = 1, k. 

Note that (ii) holds without loss generality for any bounded sequence. In- 
deed, if \a k \ < M for all i < k and k g N, then the normalized sequence (a fe )fe 
with a k = 1/2 + a k /2M satisfies (ii). 

Example 1: Let (ak)k be any convergent sequence in [0, 1]. If we let a k = a,; 
for i < k, then it is easy to see that (a fe )fc 6 S. 

Example 2: As an example of a sequence (& k )k € S that is not necessarily 
obtained from a convergent sequence, consider a bounded sequence (a fc )fc that 
satisfies for any < 8 < 1: 

#{1<J <k:\a k +1 -a k \=Q(k- s )} 

J — s-l as k — s- oo. 

k 

One can easily verify that (& k )k satisfies condition (i) of Definition 12.11 For 
instance, any uniformly distributed sequence (a fe )fe modulo one with discrepancy 



Dk = max 

i<i<fe 

falls into this category [TP] , 



k j- 1 
a.; : — 



= 0(k~ 5 ) 



Our first trace formula is concerned with the moments of Jacobi matrices 
J(a fe ,b fe ) whose defining sequences (aL k )k and (b k )k belong to S. 

Proposition 2.2. Let (a k )k, (b fc )fc be two sequences in S. For any n 6 N, we 

have 

L«/2J | k 

Trace [J"(a fe ,b fe )] = ^ { .^_ ^ }| g(of T^^f* + (fc). (2.1) 

Proof: We write Jfc := J(a fc ,b fc ) as the sum of three matrices 

J fc = L fc + D fe + 
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where Dk = diag(aj, a|) and Lk is the lower triangular matrix given by 





( 











••• \ 
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(2.2) 
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For any n e N, J£ = (L^ + Dk + Lk) n is the sum of 3™ matrix monomials 
of the form 

A 1 A 2 ---A n 

with Aj being either Lk, Dk or L\. One can easily verify that the monomials 
having a different number of Lk and L\ in their expressions have zero trace. 
Therefore, it suffices to consider those having the same number of Lk and L\. 

The assumptions on the sequences (a k )k and {h k )k imply that when we 
permute any two consecutive matrices in the product A\A 2 - ■ ■ A n , the trace of 
the resulting expression differ from the trace of the original one by o(k). In 
other words, we have 

i Traced • • • A,A l+1 ■ ■ ■ A n ] = -j-Tracc[Ai • • • A i+1 A t ■ ■ ■ A n ] + o(l). 
k k 

To see this, let us consider the case where A^ = L\ and A i+ i = L k ; the other 
cases can be handled in a similar manner. We have that 

A 1 A 2 ---A n = A X A 2 ■ ■ ■ Ai +1 Ai ■ ■ ■ A n 

+AiA 2 ■ ■ ■ A l - 1 [A l ,A i+1 ]A l+2 ■ ■ ■ A n . 

By an elementary property of the trace, we also have 

Traced A 2 ■ ■ ■ A i _ 1 [A u A l+1 ]A i+2 ■ ■ ■ A n ] 

= Tracep;, A i+1 ]A l+2 ■ ■ ■ A n A x ■ ■ ■ A^ x ] 

The product Ai+ 2 • • • A n A\ ■ ■ ■ Ai-\ is a matrix whose diagonal elements are 
bounded, while the commutator matrix [^4j,Aj+i] is a diagonal matrix whose 
elements are (b k ) 2 , (b 2 ) 2 - (b k ) 2 , (^) 2 -(&fc_i) 2 , (bl) 2 . Under the assumptions 
that {si k )k and (b fe )fe belong to S, we easily deduce 

fe-i 

Trace[[A i ,A i+1 ]A i+ . 2 ■ ■ ■ A n A x A 2 ■ ■ ■ A^} < \b k \ 2 + £ \b k +1 - b k \ 2 + \b k \ 2 

= 0(k). 

Consequently, it suffices to consider terms of the form (LkL^) 23 D^ 2 -' for which 
the trace can easily be computed. Indeed, we have 

T ra ce[(LkLl) 2 Wr 23 } = YMT~ 23 ($? J - 

i=i 
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Since there are fjrjT^^jTj monomials containing j matrices Lk and and 
n — 2j matrices Dk , we finally obtain 

i Trace Lr(a fc ,b fc )l 

K 



fc ^ (i0 2 (n-2i)! 

1 h/2j . ft 

j=0 w ' v J; i=l 

as desired. □ 



The boundedness condition (m) in Dehnition l2 . 1 l ean be removed if we slightly 
strengthen the first condition. More precisely, we have: 

Definition 2.3. We say that the sequence (a )k belongs to S' if it satisfies the 
following two conditions: For any S € (0, 1), 

(i') E?=i \a k i+ i ~ a "l = °( kl ~ 5 ) ™k^oc, 
(ii') max |ffl*| = O(logfc) as k — > oo 

l<i<k 

It is readily seen that a similar proof as the one of Proposition 12.21 holds 
when the sequences (ar)k and (b fe )fe are both in S' . We state this result as our 
second trace formula. 

Proposition 2.4. Let (ar)k, (b fe )/c G S' . For any n £ N, we have 

L"/2J , fe 

Trace [j«(a fe ), (b fc )] = £ (j , )2( ^ 2j) , + *(*)• (2-3) 



3 Main trace formulas: The //-distributed case 

We now derive our main asymptotic trace formulas for Jacobi matrices whose 
defining sequences (a and (b fc )fc are distributed according to some probability 
measure. We consider different cases depending on if the sequences (a fc )fc and 
(hk) k are bounded, unbounded or monotone. 



3.1 The bounded case 

Our first results are concerned with compactly supported probability measure 
/j on I 2 . After normalization, we may assume without loss of generality that 
the support is contained in I 2 :— [0, 1] x [0, 1]. 
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Definition 3.1. We say that a sequence (a k ,h k )k with a k and b k in [0,1] for 
all j = 1, k is /i- distributed if for any continuous function ip on [0, l] 2 , one 
has 



, lim l5Z^( a j' 6 j) = / 4>{x,y) dn{x,y). 

fc-S-OO k J J Jj2 

Example 3: Consider the sequences (a fc )fe and (b k )k that are convergent, i.e. 
a fe = (ai,...,afc) and b fc = (bi,...,bk) with au — > a and bk — > b. It is easy to 
see that (a k ,b k )k is /i-distributed with \i — S a x 8b, the product of the Dirac 
measures at a and at b. 

Example 4: Let a, b : [0,1] — > [0,1] be two continuous functions. We define 
the sequences a k and h k by 

a fc = (a(l/k),a(2/k), ...,a((k - l)/k),a(l)) 

and 

b fe = (b(l/k), b(2/k), b((k - l)/k), 6(1)). 

It follows that (a k ,b k ) k is to /-distributed where to/ is the probability distribu- 
tion associated to the random variable f{x) = (a(x), b(x)), i.e. 

tp(x,y) drrif(x,y) = / ij)(a(x), b(x)) dx 
[o,i] 2 Jo ' 

for any ip G C{I 2 ). 

In our first theorem of this section, we give a trace formula for Jacobi ma- 
trices with bounded sequences that are /^-distributed. 

Theorem 3.2. Let [i be a probability measure on I 2 . Suppose that (a k ,h k )k is 
a /j,- distributed sequence with (& k )k, (b fc )fc G S. Then, we have 



\imyTrace[c/)(J(a k ,b k ))} = -[ [ (f>(x + 2y cos t) dt 

k^-oo k TT Jj2 J 



dfi(x,y) (3.1) 



for any <f> G C[— 3, 3] . 

Proof: For any nonnegative integer n, Proposition 12.21 implies 

\ Trace [j"(a fc ,b fc )] 
k 

■E^/^V*,),.?). (3-2, 

In order to replace the sum X)}=o by the sum over all multi-indices a = 
(ol\, ct2, 0:3) with |a| = n, we introduce the function sine : R — ¥ R defined by 

sinc(0 = ^ J e***dt. 
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Note that sinc(O) — 1 while sinc(7ra;) = when X is a non-zero integer. Using 
this function, (|3.2I) becomes 



1 



Trace [j"(a fe ,b* 



If f E (^)e™ t(ai - a2) x a3 y ai+a2 dfi(x,y)dt + o(l). 



\a\—n 

The multinomial theorem and a simple change of variables then yield 
1. 



; Trace [J"(a fe , b fc )] 



i 2 Jo 



(x + 2ycos(t)) n dtdfi(x,y) + o(l). (3.3) 



This establishes the result for the moments of J(a fc ,b fc ). By linearity, the 
result also holds for polynomials of arbitrary degree. Now let <j> be a continuous 
function on [—3,3]. Note that 

{x + 2ycos{t) : (x,y) <= I 2 , t £ [0,vr]} C [-1,3] 

and cr(J(a fc ,b fc )) C [— 2\/3, 2\/3]- By Weierstrass Approximation Theorem, 
there is a polynomial P such that \\<f> — P\\oo < e /3 on [—3, 3]. In particular, it 
implies that 



iTrace[0(J(a fc ,b' £ ))] - -Trace[P(J(a fe , b fe ))] 



<e/3, 



(3.4) 



together with 



j 2 Jo 



0(a; + 2y cos(t)) dt dfj,(x, y) 

P(x + 2y cos(t)) dt dfj,(x, y) 



i 2 Jo 



< e/3. (3.5) 



Finally, (|3.3|l shows that we can choose k large enough so that 



\^[P{J^))]-\f f p 



(x + 2ycos(t)) dtdjjb{x,y) 



< e/3. (3.6) 



The conclusion of the theorem follows by combining (|3.4I) , (|3.5p , and (j3.6|) . □ 

Example 5: Let (a fe )^ and (b fc )i, be convergent sequences as in Example 3. 
From Theorem I3.2[ we easily deduce that 



lim -J-Tracc[0(J(a fc ,b fe ))] = - f 4>(a + 2b cos t) dt 

k— >oo k 7T Jo 



a+2b 



n Ja-2b 



dt 



u — a + 2b) 
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This type of Jacobi matrices naturally arise for orthogonal polynomials that 
belong to the M(a,b) class introduced by Nevai [11]. One can use this result 
to easily compute the asymptotic distribution of the zeros of many classical 
orthogonal polynomials such as the Jacobi (a = 1/4, 6 = 1/2), Chebyshev 
( a = 0, b = 1/2), Legendre (a = 0,6 = 1/2) and Gegenbauer (a = 0,6 = 1/2) 
polynomials. 

Example 6: In computing the asymptotic distribution of the zeros of Van 
Vleck polynomials [4], one is led to consider a Jacobi matrix of the form 





( 


a i 


6i 








\ 








6i 


«2 


6 2 








Jk — 







l>2 


A3 







(3.7) 












bk-i 








\ 










6fc-i afc 


/ 





where the entries are up to some constants given by 

(j - 1) 
a,; = ; and 6, 



It is easy to see that the sequences (a,k)k and (6fc)fc are both in 5, and the 
sequence (a,k,bk)k is /i-distributed with 



(j-{x,y) 

It then follows from Theorem 13.21 that 



5 (y - - x)) 



2v^ 



1 1 f 1 f" 

lim — Trace [(/>(Jfc)] = — / / (j)(x + 2\/x(l — x) cost) dt 

k->oo k 7T J J 



dx 



In a recent paper, Kuiljaars and Serra-Cappizano [8 proved a general result 
regarding the asymptotic distribution of the eigenvalues of Jacobi matrices. 
Their results generalize earlier ones obtained by Kuiljaars and Van Assche [9] 
and Geronimo, Harrell II and Van Assche [7]. We give here a new proof of 
their result based on our previous trace formula in the generic case \a!* \ < 1 and 
< < 1. 

Corollary 3.3. (Kuiljaars- Serra Cappizano) Let (J(a fc ,b fe ))j, be a sequence of 
Jacobi matrices and suppose we can choose two bounded measurable functions 
a, b : [0, 1] — > [0, 1] satisfying the conditions: For any e > 0, 

|{s e [0, 1] : \a\ sk] - a(s)\ > e}\ ->• as k -> oo, (3.8) 
\{s e [0, 1] : |6* - 6(s)| > e}| -J- as k -> oo. (3.9) 
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Then, we have 

i i f 1 

lim - Trace[(j){J{a k ,h k ))] = - / <f>{a(x) + 2b(x) cosy) dxdy 

fc->oo k 7T J J Q 

for any tf> 6 [—3, 3]. 

Proof: For k large enough, conditions (|3.8I) and (|3.9[) together with Lusin's 
Theorem imply that we can choose the functions a and b to be continuous on 
[0, 1] except for a set of measure 1/fc and for which the additional conditions 

#{j<k:\a*j-a(j/k)\>e} = o(k) 

#{j<k:\b k j -b{j/k)\>e} = o{k) 

hold for any given e > 0. It easily follows that for any continuous function ip on 
[0,1] 2 , 

lim Jy^( ffl k ,b k ) = { Ma(x),b(x))dx. (3.10) 

If we introduce the map / : [0, 1] — > [0, l] 2 defined by f(x) = (a(x), b(x)), then 
it follows as in Example 4 that the sequence (& k ,h k )k is m^-distributed. The 
conclusion is then an immediate consequence of Theorem 13.21 □ 



In the next result, we derive a simple extension of Example 5. Indeed, we 
consider sequences (si k )k and (h k )k that have a finite number of accumulation 
points. 

Corollary 3.4. Let (a k )k and (h k )k be two sequences in S . Let SI = \uj\, ...,o; n ] 
be an n vector for which < Wj < 1 and 



\n\ :=^^ = l. 



Suppose there exists two n-tuple ...,a n ) and (j3±, /3 n ) such that for any 
1 < p < n, and e > 0, 

^■■\a k ~a P \<eand\b k -b p \<e} ^ ^ ^ 

Then, for every continuous function <f> on [—3,3], one has 

k 

[x + 2y cost) du>(x, y) dt (3.12) 



lim i V Tm Ce [0(J(a fc ,b fe ))] = - f / 



where u> is the probability measure given by the weighted sum of delta measures, 
i.e. 

n 

lo(x, y) = y^ j u i 6 ai (x) x 6 bi (y). (3.13) 
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Proof: It is readily seen that both sequences satisfy the assumptions of Theorem 
I3.2l with fi given by uj as defined in the statement of the theorem. □ 

As a simple consequence of the previous result, we have the following result 
also due to Kuiljaars and Serra Capizzano [8] when the sequences (a k )k and 
(b fc )fe have a single accumulation point. 

Corollary 3.5. (Kuiljaars-Serra Capizzano) Let (<Zfc)fe and (6fe)fe, bk > 0, be 
two bounded sequences as in Example 1, i.e. (sL k )k is the sequence defined by 
a fe = (ai, ak) and similarly for (h k )k- Suppose there exist real constant a and 
b > such that for every e > 0, 

#{j < k : |aj - a\ > e} = o(k) as k -> oo, (3.14) 

and 

#{j < k : \bj - b\ > e} = o(k) as k -> oo. (3.15) 
Then, for every continuous function <f> on [a — 26, a + 2b], 

lim i V Trace[0(J(a fe ,b fe ))] = i / 0(a + 2&cosx) tte. (3.16) 

Remark: Recently, Trench [12] obtained weaker conditions under which the 
conclusion of Corollarv l3.5l holds. However, his conditions are of different nature 
as they are related to the spectrum of J(ak, bk) rather than the sequences (afc)fe 
and {b k )k- 

3.2 The unbounded case: Part I 

We now turn our attention to the unbounded case, i.e. we now assume that the 
sequences (a fe )^ and ih k )k can take values over the whole real line. Our first 
result is concerned with sequences that can be contracted in order to fit the 
bounded case. The next definition is based on the terminology introduced by 
Van Assche [13] . 

Definition 3.6. We say that r : (0, oo) — > (0, oo) is a regularly varying function 
(at infinity) for the sequence (a k , h k )k if the normalized sequences 

( a fe \ / b fc \ 

are in S . 



The next result is an immediate consequence of Theorem 13.21 so its proof is 
omitted. 



11 



Corollary 3.7. Let r be a regularly varying function for the sequence (a ,b )fc. 
Suppose that and (a k /r(k), h k /r(k))k is \i- distributed for some probability mea- 
sure fi on I 2 . Then,, we have 



lim — Trace 

k-t-oo k 



for any <f> <G C[— 3, 3] 



r(k) ' r(k) 



J I 2 JO 



x + 2y cos t) dt dfi(x, y) 



Example 7: Using the three-terms recurrence relation for the Hermite poly- 
nomials, it is well-known that the zeros z k ,...,z% of the fcth degree Hermitt 
polynomial are the eigenvalues of the tridiagonal matrix 

/ 1/2 ••• \ 
1 1/2 '•■ : 
2 '•• 

: ■■. ■■. ■■. 1/2 

V ••• k — 1 / 
It easy to see that Hk is similar to the Jacobi matrix 
/ 1/V2 



H k = 



(3.17) 



o \ 



1/V2 




V 



(3.18) 



Indeed, H' k = DkHkD^ where Dk = diag(di, dk) is the diagonal matrix 
whose elements are defined recursively by the relations 

d\ = 1, and d i+ i = -^= for i = 1, k — 1. 
\/2i 



The regularly varying function is r{k) = \2k and the sequence 

k 



i k b k 



r(k) 1 r(k)J k 







k 



is ^-distributed on I 2 with 



dfi(x, y) = 2yS (x) dxdy. 
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Consequently, the asymptotic distribution of the contracted zeros of Hermite 
polynomials is given by 



lim > i 

■—v™ ^ — < 



1 



2k 



n Ji 2 Jo 
2 

7T 



JO 
2 



c/>(x + 2y cos t)dt dfi(x, y) 

4>{2y cos t) y dt dy 

= — I <t>{x) \J A - x 2 dx 
2tt J_ 2 

for any <f> G C[— \[2, V2]. We thus obtain another proof of the well-known fact 
that the asymptotic distribution of the contracted zeros of Hermite polynomials 
is the Wigner semicircle distribution. 

Example 8: One can use a similar approach to derive the well-known fact that 
the asymptotic distribution of the contracted zeros of Laguerre polynomials is 
the Marchcnko-Pastur distribution. Indeed, the zeros z\,...,z\ of the fcth degree 
Laguerre polynomial are the eigenvalues of the Jacobi matrix 



Lk 



( 1 
1 



V o 







fc-1 

fc-1 2k -\ J 



(3.19) 



Hence, the regularly varying function is r(fc) — k and the sequence is 



2,7 - 1 
k 



3=1 



3 = 1/ 



j(k)' r(k)J k 

is ^-distributed on I 2 with 

d^(x, y) =S {x- 2y) dxdy. 



The asymptotic distribution of the contracted zeros of Laguerre polynomials is 
therefore given by 



lim > < 

k— ¥oo — ' 



3 = 1 



1 

I" J J- ./() 
1 

T JO JO 



4>{x + 2y cos t) dt dfi(x, y) 
<f>(2y + 2y cos t) dt dy 



dx 
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for any <p g C[0,4]. 

More generally, we can use our trace formula in Corollary 13.71 to give a new 
proof of the following result due to Van Assche. 

Theorem 3.8. Let r be a regularly varying function with parameter a € (0, 1), 
i.e. for every t > 0, 

x^co r(x) 

Let (ak)k o,nd {bk)k be two sequences with ak S M and bk > satisfying the 
conditions 

lim — = b. 



lim -77T = a, 



fc^oo r(fc) 



fc-s-oo r(k) 



(3.20) 



If J(& , b ) is £/ie Jacobi matrix with a = (oi, as) and b = (&i, i/ien 
we /iaue 



lim — Trace 

k— s-oo k 



Si h k 



r(k) r(k) 



1 



a+2b 



f(x) v a (x)dx. 



a-2b 



Here v a (x) denotes the Nevai-UUman density with parameter a defined as the 
Mellin convolution 

f 1 dy 

v a (x) = b a * uj a ,b{x) = I b a (y)uj a ^(x/y) — 

Jo v 



with 



b Q (y) 



ay"- 1 if ye (0,1), 







elsewhere, 



J afi (y) = I W ^(a+2b-y)(y~a+2b) 

I elsewhere 



if y^ (a-2b,a + 2b). 



Proof: The conditions (|3.20p imposed on the sequence (ak)k and (bk)k imply 
that for any e > 0, 



#{l<J<fc: 



r(t) 



)°l<'-} 



and 



#{l<j<* = \^~b(j:) a \<e} 



for k — > oo. In particular, it follows that 



lim —Trace 

k— ¥oo k 



r(k) ' r(k) 



= lim -Trace U(j(au fe , &u fc ))l 
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where the sequence u fc is given by 

u fc = (l/fc Q ,(2/fc)« ! ...,((fc-l)/A : )«,l). 

Clearly, u fe £ S and by Example 4, the sequence (au k , &u fc ) is /i-distributed 
with fi satisfying 

/ ip(x, y) dfx{x, y) = / ijj(ax a ,bx a ) dx 
Ji 2 Jo 

for any ip £ C(I 2 ). By Theorem 13.21 we deduce that 



lim —Trace 

fc->00 k 



r{k) ' r(k) 
1 f 1 



o Jo 



(j)(ax a + 2bx a cos(?rf)) dt dx 



1 f 1 f a+2b „ , dxdt 



f(t a x 



o J a -2b \/{a + 2b — x)(x — a + 2b) 

Under some simple changes of variable, last is easily seen to be equivalent to 
the expression given in the statement of theorem. □ 

3.3 The unbounded case: Part II 

Our second result for Jacobi matrices for unbounded sequences is concerned with 
sequences (& k )k and (h k )k that belong to S' (see Definition [23]) . Before we can 
go further, we need to reformulate Definition 13 . 1 1 for unbounded sequences. 

Definition 3.9. We say that a sequence (a k ,h k )k with (a k )k and ih k )k in S' 
is fi- distributed if for any if) £ Cb(M 2 ) := C(R 2 ) n L 00 (IR 2 ) ; the space of bounded 
continuous functions on R 2 , one has 

lim \ y^^{a k ,b k )= f ip(x,y) dft(x,y) 

for some probability measure pi onM 2 . 

Recall, a probability measure is said to be tight if for any e > 0, one can 
find a closed bounded rectangle TZ = x [0:2, A2] such, that 

A*(M 2 - TV) < e. 

Theorem 3.10. Let (a fe ) fe and (h k )k be two sequences in S' . If (a. k ,h k )k is 
pi- distributed for some tight probability measure fx on M 2 , then we have 

lim \ Trace U(J(a fe ,b fe ))l =- / / Mx + 2y coat) dfj,(x, y) (3.21) 

fe^oo k 7T J R 2 J 

for any <f> £ 
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Proof: By Proposition 12.41 we know that 

i T race[j"(a<%b*)] = - £ Qjpfi^fi E «T*ff + oft)- (3-22) 

Since /x is tight and (a fc , b fc )/j is /i-distributed, there exists for any e > a closed 
bounded rectangle RCI 2 such that 

V j; - < e 3.23 

for /c large enough. We now introduce the smooth cutoff function 9 : R 2 — > [0, 1] 
defined by 

\ i£(x,y)eK 
i{(x,y)<£K ( 



9(x,y) 



where 7Z e is a compact subset of R 2 with 1Z C lZ e and m(lZ e /lZ) < e. Equations 
fl3~22]) and (pHS]) yield 

i Trace [J"(a fc ,b fc )] 

, L«/2J , k 

= ~ k Y, {j])2{ 2'_ m B^OT^ 6.) + "(1) (3.24) 
for any nonnegative integer n. Moreover, (a fe ,b fe )fc being ^-distributed implies 
i Trace Lr(a fc ,b fc )l 

K 

L"/2J , 

= E ( J - ! )2 ( r'-2j)! J^- 2 ^ 0{x,y)diM(x,y)+o(l). (3.25) 

We now perform our usual trick, that is we introduce the sine function in (|3.25[) 
and use the multinomial theorem to obtain 

r Trace \j n (a k ,b k )] 
k 



(x + 2y cost)™ 6(x,y)dfj,(x,y) + o(l). (3.26) 



Note that 

K\ := {x + 2ycost : (x, y) £ lZ e and t € [0, 7r]} 

is a compact subset of R. Moreover, the assumptions on fi together with (|3.23[) 
imply that there exists a compact subset C R such that 

#{A : A e a( J(a fc , b fc )), A £ tf 2 } = o(fc) 



1G 



as gets k arbitrary large. In particular, we have 



iTrace[</"(a fc ,b fc )] =\ ]T A™ + o(l) 



(3.27) 



\£a(J(a k ,b k ))DK 2 



By Weierstrass Approximation Theorem applied tp a closed bounded interval 
that contains K\ U K% , it is then easy to see that (|3.26[) yields 



1 



-Trace [^(J(a fc , b k ))] 



1 



(x + 2y cos t) 9(x, y) d/i(x, y) + o(l) (3.28) 



n Js. 2 Jo 

for any <f) G Cfe(K). Furthermore, since the measure /x is tight, it follows that 
i>(x + 2y cost) (1 — 6(x, y)) d/i(x, y) 



< 



>R € -R 

< 2U\\ no e. 



(x + 2y cost) dfj,(x, y) 



(x + 2y cost) (l-0(x,y))dfi(x,y) 



(3.29) 



The desired trace formula (13.21| is then an immediate consequence of (|3.26j) 
and ff3T29|) . □ 

Example 9: Let (X fe )j. and (Y k )k be two sequences with X fc and Y fe being 
fc-vectors whose components are i.i.d. random variables that are normally dis- 
tributed with mean zero and variance a 2 (k) — > as fc — >oo. First, we note that 
for all l<j<k, 



Probjw E K : \X k (uj)\ > logfc} 

1 / -log 2 fc/2<x 2 (fc)' 

= ^ Er fc(.( fc) log fe )= ^ fclQgfc 

Secondly, Chebyshev's inequality implies that 

Prob {wet: \X!? +1 (u) - Xf(u)\ > <r 1/2 (fc)} 



(3.30) 



< 



1 



< 



a(k) 
8 

8a(k) 



2tt 



-A, ,\\2 -uj 2 /2a 2 (k) 



2n 
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From this, we easily deduce that 

k 

i £ - *»l = cV' 2 (fc)) (3-31) 

3=1 

for asymptotically almost every (a.a.e.) to G M. In particular, if we assume 
that a(k) = 0(k~ s ) for some 5 > 0, it then follows from (pHUl) and (pOTj) that 
(X fe (u;))fc € 5' for a.a.e. w € R. Of course, a similar conclusion also hold for 

(Y fe M) fc . 

Furthermore, under the independence assumption, one can apply the strong 
Law of Large Numbers to conclude 

k 

i ^-M) = 2^(fc)l 2 V ' (x ' 2/) ^ (a;2+? ' 2)/2CT2(fc)dxd2/ + o(1) 

= V(0,0) + o(l) 

for a.a.e. u£R and any -0 6 Cfc(R 2 ). In other words, the sequence (X k (uj), Yfc(w))fc 
is <5o x ^-distributed in the sense of Definition 13.91 By Theorem 13.101 we con- 
clude that for a.a.e. wet, 

lim -Ttace[0(J(X fc (w),Y*(w)))] =0(0) 

for any G C 6 (M). 



3.4 The monotone case 

An obvious, but interesting corollary of Propositions 12.21 and 12.41 is obtained by 
considering sequences that are almost everywhere monotone. 

Definition 3.11. We say that (a fc ) fceN is a.e. increasing if 
#{l<i<k:a k < a* } > 1 



k 



as k — ¥ co. 



Of course, a similar condition holds for a.e. decreasing sequences. Clearly, 
any a.e. monotone sequence that satisfies conditions (ii) or (ii') of Definitions I2.1 1 
and 12 .3! respectively will also satisfy conditions (i) or (i') of the same definitions. 
As a consequence, we obtain the following result. 

Proposition 3.12. Let (a fc )fc, (h k ) k be two a.e. monotone sequences with either 

(a k ) k , (b k ) k e S 



or 



( a k ) k ,(b k ) k eS'. 
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Then, for any neN, we have 

L«/2J 

Trace [,J n (a k ,b k )] = ^ 

3=0 



{j\f{n-2j)\ 



(a k r-^(b^+o(k). (3.32) 



If we assume furthermore that the sequence (a fc , h k )k is ^-distributed as in 
the sense of Definitions 13. ll or 13.91 then we easily derive the following two results 
that are simple consequence of Theorems 13.21 and 13.101 

Corollary 3.13. Let (a fc ) fc and (h k )k be two a.e. monotone sequences in [0, 1]. 
// (a fe , h k )k is [i- distributed for some probability measure ji on I 2 , then we have 



lim i Trace U(J(a fc , b fe ))l =- f f , 

fc->oo k J 7T Jj2 J 



(x + 2y cost) dfj,(x, y) (3.33) 



for any (j) 6 C[— 3, 3] . 

Corollary 3.14. Let (a fe )fc and (h k )k be two a.e. monotone sequences in K 
that satisfy the estimates 

\a k \ = O(logfc) and \b k \ = O(logfc) for all 1 < j < k, k e N. 

If (& k ,b fe )fe is ^l- distributed for some tight probability measure pi on M. 2 , then we 
have 

lim \ Trace U(J(a fe , b fe ))l =- / / <p(x + 2y cost) dfi(x, y) (3.34) 

k^oo k 7T J R 2 J 

for any <f> 6 



Example 10: Let (X fe ) fc and (Y k ) k be sequences where 

X k = (X k , ...,X k ) and Y fc - (Y k , ...,Y k ) 

are two k random vectors. We assume for all 1 < j < k, k € N, that the X k, s 
and Y k, s are i.i.d. random variables of mean zero and variance one on some 
probability space (fi, P). We denote by X*- fe ' the random vector whose jth entry 
is the jth order statistic Xk-. obtained from X k , ...,X k , i.e. 

X^=(X k 1} ,...,X k k) ) 

and similarly for Y^ k \ For any w £ fi, we denote by Jfe(w) the Jacobi matrix 
given by 

( x k {1) {u) y ( fe 1} H o ••• o \ 

YftM X( 2) { U ) Y k 2) {u) \ 
J(fe)M= Y k 2) {u) X k 3) (uj) ■■. • ( 3 - 35 ) 

'■• '•• '■• 
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By definition of order statistics, the sequences X( fc ) and Y^ k ' are monotone 
increasing. If we assume that X k (oj) G [0, 1] and Y k (uj) G (0, 1], then the strong 
Law of Large Numbers implies that for a.a.e uj G f2, 

lim 5^V(-^y)H.ly)M) = / tfav) dfj,x,r(x,y) 

3=1 7 

for any ip G C(I 2 ). Here, we denote by /ix,y the joint probability distribution 
of X and Y where X and Y are two random variables on fi that follow the same 
distribution as X k and Y k respectively. Consequently, we can apply Corollary 
13.131 to conclude that 

i i r r 

lim — Trace [^(Jfc(o;)] = — / / (f>(x + 2y cost) d^ix y(x, y) (3.36) 

fc->oo k TT J j2 J 

for any <f> G C[— 3, 3] and a.a.e w G il. 

Similarly, if we assume that X k — O(log k) and K* — O(log k), then we can 
apply Corollary 13. 141 to obtain 

i i r r 

lim — Trace [0(Ji;(aj)] = — / / <j>(x + 2y coat) dfix y {%, y) (3.37) 

fc->oo k TT J R 2 J 

for any </> G Cf,(K) and a.a.e w G f2. 



4 Existence Results 

In this section, we would like to address the following question: Given two 
sequences (a. k )k and (h k )k in S or S' , does there exist a probability measure 
with support in I 2 or a tight measure probability measure /x with support 
in R 2 for which the sequence (a fc ,b fc )fe is /^-distributed? It is easy to construct 
examples of sequences that are not /i-distributed for any probability measure. 
For instance, take the sequence (a k ,h k )k with 

r „fc^ _J(0,1) if2|fc 
if 2 f fe. 

However, as the next result shows one can always find subsequences that are 
/^-distributed. 

Proposition 4.1. Let (& k )k and (b fc ) fc be two sequences in S. Then, there 
exists a subsequence of positive integers kj and a probability measure fi with 
support in [0,1] 2 such that (a 3 ,b ')j is /i- distributed. 

Proof: For any x, y G ffi, consider the distribution functions Fk defined by 

#{j < k : a) < x and b k < y} 
Fk{x,y) = J — * . 
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By Helly's selection principle ([3], Theorem 25.9), we can find a subsequence of 
positive integers kj, kj — > oo as j — > oo, such that Fk, converges weakly (in the 
probability sense) to a distribution function F. By Theorem 12.5 in [3], we can 
associate to F k . and F unique probability measures /i^ and \x on [0, l] 2 defined 

by 

H((a, b) x (c, g0) = F(6, d) - F{b, c) - F(a, d) + F(a, c) 

and similarly for /i}-. . Since fx^ . converges weakly to n, we can apply Theorem 
29.1 in [3J to conclude that (& kj , h kj ) is /^-distributed as in the sense of Definition 

ED □ 

Proposition 4.2. Let (a fc ) fc and (b fe )fc be two sequences in S' . Suppose that 
for any e > 0, there exists a bounded rectangle R such that 

#{l<j<k:\tf\tR} m<J<k:\b k \^R} 

< e, and < e (4.1) 

k k 

for k large enough. Then, there exists a subsequence of positive integers kj and 
a tight probability measure on R 2 such that (a k \h kj )j is [i- distributed. 

Proof: The argument is identical as for the proof of Proposition 14.11 with the 
exception that Theorem 29.1 in [3 must be replaced by Theorem 29.3. □ 

By combining Proposition 14.11 and Theorem I3.2[ or Proposition ^. 21 together 
with Theorem I3.10[ we obtain the following two existence results. 

Corollary 4.3. Let (a fc ) fe and (h k )k be two sequences in S. Then, there exists 
a subsequence of positive integers kj and a probability measure fi with support 
in I 2 such that 

]im^-Trace[(j)(J(a k i,b k '))] = -[ [ <j>(x + 2y cos t) dt dfx(x, y) (4.2) 
J-><» kj it Jj2 J Q 

for any <f> 6 C[— 3, 3] . 

Corollary 4.4. Let (a k ) k , (b k ) k e S' be two sequences satisfying (|4.ip . Then, 
there exists a subsequence of positive integers kj and a probability measure pi 
such that 

]im^-Trace[<f>(J(a k i,b k J))] = -[ [ (j>(x + 2ycosi) dt dfx{x,y) (4.3) 

j->oo kj 7T J R 2 Jq 

for any tf> € C b (R). 

5 Conclusion 

Our results establish new and very general conditions for which one can compute 
explicitly the asymptotic distribution of the spectrum of Jacobi matrices. Our 
results are new in the sense that we do not assume any type of convergence 
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for the sequences (a k )k, (b k )k of J(a fe ,b fc ), but we only assume they satisfy a 
small-deviation condition and they are /i-distributed. 

A quick a look at the proofs should convince the readers that similar tech- 
niques as the ones developed in this paper can be used to compute an asymptotic 
trace formula for band matrices. For instance, it is straightforward to extend 
the proof of Propositions 12.21 and 12.41 to matrices with three diagonal bands of 
the form 

(4 b\ \ 

c\ a k 

a k 

V C t- P a k j 

where the sequences (a fe )fc, (b fc )fc, (c k )k are assumed to satisfy the usual small- 
deviation conditions. In particular, we would like to use these extensions to de- 
rive new results for the asymptotic distribution of difference and Sturm-Liouville 
operators. We hope to address these points in future research work. 



U k-q 



(5.1) 
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